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We consider a Lebwohl-Lasher model of chiral particles confined in a planar cell (slit pore) with 
different boundary conditions, and solve it using mean-field theory. The phase behaviour of the 
system with respect to temperature and pore width is studied. Two phenomena are observed: (i) an 
isotropic-cholesteric transition which exhibits an oscillatory structure with respect to pore width, 
and (ii) an infinite set of winding transitions caused by commensuration effects between cholesteric 
pitch and pore width. The latter transitions have been predicted and analysed by other authors for 
cholesterics confined in a fixed pore and subject to an external field promoting the uniaxial nematic 
phase; here we induce winding transitions solely from geometry by changing the pore width at zero 
external field (a setup recently explored in Atomic-Force Microscopy experiments). In contrast with 
previous studies, we obtain the phase diagram in the temperature vs pore width plane, including the 
isotropic-cholesteric transition, the winding transitions and their complex relationship. In particu¬ 
lar, the structure of winding transitions terminates at the capillary isotropic-cholesteric transition 
via triple points where the confined isotropic phase coexists with two cholesterics with different 
helix indices. For symmetric and asymmetric monostable plate anchorings the phase diagram are 
qualitatively similar. 


I. INTRODUCTION 

Cholesteric phases of liquid crystal materials are im¬ 
portant from fundamental and technological points of 
view 0,0, and a lot of effort is still being expended 
in their study. Interest in cholesterics stems from the 
existence of a direction, perpendicular to the nematic di¬ 
rector, about which the latter winds up with a specific 
periodic length, the pitch, causing the material to be op¬ 
tically active at optical wavelengths. The pitch sensitiv¬ 
ity to temperature and external fields, and the ensuing 
varying reflectivity properties, can be exploited in appli¬ 
cations. Also, temperature variations of the pitch can be 
suppressed by different techniques, e.g. by mixing two 
chiral dopants with opposite temperature dependence in 
their pitch 0. 

Cholesteric materials continue to be of interest, in par¬ 
ticular, as the basis for the design of the twisted nematic 
cells of liquid-crystal-based devices. Many experimental 
and theoretical works have been devoted to this problem. 
Still many basic issues remain unexplored in the field of 
capillary phenomena in cholesteric materials. For exam- 
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pie, the first theoretical study of the isotropic-cholesteric 
interface is very recent 0. 

Most theoretical studies of cholesterics have been based 
on the Landau-de Gennes model, especially in connec¬ 
tion with the existence of the blue phases (see e.g. i). 
Another fruitful approach is the chiral Lebwohl-Lasher 
model @0, which has been studied from various per¬ 
spectives. The bulk properties of the chiral Lebwohl- 
Lasher model were explored by van der Meer et al. 0 
using mean-field theory. The most striking conclusion 
was that the pitch was independent of temperature, and 
that higher-order terms in the interactions were necessary 
to restore a temperature-dependent pitch. The model 
was first simulated by Saha et al. [3 using the Monte 
Carlo technique and a planar rotor model for the par¬ 
ticles. Discrepancies were found initially concerning the 
independency of the pitch with temperature. Later sim¬ 
ulations by Saha and Saha 0 of the full model (without 
orientational restrictions) and of the planar rotor model 
by Luckhurst et al. 0 confirmed the mean-field predic¬ 
tions. The latter authors obtained a continuous isotropic- 
cholesteric transition, and explained the temperature- 
independent feature as a result of the existence of a he¬ 
lical transformation, suggested by T. J. Sluckin, that re¬ 
moves the chiral term in the energy. Note that, in more 
realistic models with density fluctuations, one expects the 
pitch to depend on the temperature due to variations in 
density and order parameter with temperature. 
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Other authors have used the model to explore dif¬ 
ferent features of chiral phases. Memmer and Janssen 
El studied the influence of the fourth-order terms on 
the pitch, using simulation and self-determined bound¬ 
ary conditions. In agreement with van der Meer et al. 
0, temperature-dependency of the pitch was restored. 
A helix inversion at a specific temperature was obtained 
by these authors with the same fourth-order interaction 
potential [l^ , again in accord with the mean-field calcu¬ 
lation of van der Meer et al. Memmer and Janssen 
El also studied the induction of a helical twist by mix¬ 
ing chiral dopant and a guest achiral system, a common 
procedure in experiments and applications to optimise 
device performance. 

The chiral Lebwohl-Lasher model has also been inves¬ 
tigated in confined geometry. This is a relevant model 
to study twisted liquid-crystal displays. These systems 
operate with the help of an external field. To cite a 
representative example of a simulation analysis, Mem¬ 
mer and Fliegans studied a chiral Lebwohl-Lasher 
model confined by two parallel plates and subject to 
a constant external field normal to the plates. A cor¬ 
responding Freedericksz transition between twisted and 
uniform director arrangements was characterised. An ex¬ 
ternal electric field can be used to control the pitch in 
a bulk cholesteric sample. In the case of a cholesteric 
confined between two plates a fixed distance apart (pore 
width), there is a commensuration issue, and the vari¬ 
ation of the number of pitch lengths can only occur in 
discrete jumps involving structures with different num¬ 
ber of helix indices (number of half turns of the director 
about the helix direction) Ell- These variations can be 
induced by temperature or by an electric field. Previous 
and more recent works have analysed the problem using 
simple elastic models [TtI - I^ . There are various experi¬ 
mental works that deal with these problems (25l - [3^ . 

Our purpose in the present work is to analyse commen¬ 
suration effects in confined chiral liquid crystals when 
the agent that changes the number of pitch lengths is 
not temperature or an external field, but the pore width. 
Therefore, we are inducing the winding transitions by 
geometry alone. An important feature of our theoreti¬ 
cal technique is that order parameters and a free-energy 
functional are used; therefore the nematic and cholesteric 
order are analysed consistently as a function of tem¬ 
perature, and the confined order-disorder transition and 
structural phenomena in the slit pore can be followed 
and described collectively. Also, we consider a free- 
orientation model, so that the main axes of particles are 
not restricted to lie on planes parallel to the confining 
plates. 

To help focus our work, we can think of a material con¬ 
sisting of a mixture of two chiral substances such that the 
pitch temperature dependence is suppressed Q. This 
material can be placed in a cell with strong anchoring 
conditions, such that the pore width can be controlled, 
as in an Atomic-Force Microscopy (AFM) experiment; 
indeed such a setup has been used to identify capillary 




FIG. 1. Setup for the Lebwohl-Lasher cell: a simple-cubic 
lattice sandwiched between two parallel plates normal to the 
[100] direction. Planes parallel to the plates are indicated by 
an index from 1 to N. The pore width is h = {N -|- l)a, 
where a is the lattice parameter. The favoured alignment at 
each plate is given by unit vectors n and m, respectively, 
with a the angle between the two vectors. All lattice sites are 
occupied by one particle with orientation d>. Cartesian axes 
used (lab reference frame) are shown at left, along with the 
polar angles associated with a> (note that the z axis lies in 
the plane of the plates whereas the y axis is normal to the 
plates). 


nematization transitions and, more recently, even to 
obtain the pitch in cholesteric materials with high accu¬ 
racy [^ . The control parameter is h/p, where h is the 
pore width and p the cholesteric pitch; this parameter 
dictates the number of cholesteric periods or helix index 
in the cell, and its variation will induce transitions be¬ 
tween cholesteric slabs of different numbers of half-turns 
of the director (winding transitions). As in the AFM 
experiments (^ . an inhnite collection of winding tran¬ 
sitions are obtained. Using the chiral Lebwohl-Lasher 
model, we analyse these transitions in detail, locate the 
capillary isotropic-cholesteric phase transition as a func¬ 
tion of the pore width, and establish the connection be¬ 
tween this transition and the winding transitions, a re¬ 
lationship that previous studies have not addressed. An 
interesting topology in the temperature vs. pore width 
phase diagram is observed. The cases of monostable sym¬ 
metric and asymmetric surface alignmentsare considered, 
and the differences between the two cases discussed. 


II. THEORY 

The model used to explore the conhned cholesteric ma¬ 
terial is the chiral Lebwohl-Lasher model. This model has 
some advantages over the popular Landau-de Gennes or 
elastic models used in previous studies. The most im¬ 
portant of them is the fact that, in contrast to the latter 
approaches, the model has a microscopic nature and the 
interaction potential appears explicitly. Different terms 
in the interaction potential can be tuned to model spe¬ 
cific features of real molecules. By contrast, the Landau- 
de Gennes or elastic theories are meso- and macroscopic 
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in nature, aim at a description at larger length scales, 
and become the logical choices when the nematic direc¬ 
tor is spatially deformed on a scale much larger than 
the molecular scale, or surfaces induce the presence of 
defects. These scales are difficult to access using a mi¬ 
croscopic approach due to computational limitations. 

In the Lebwohl-Lasher model, each site of a simple- 
cubic lattice of lattice parameter a is filled with a uniaxial 
particle described by an orientation d> = {0, ip). Particles 
interact via a nearest-neighbour potential that depends 
on the two particle orientations a>, and possibly also 
on the relative position vector r. In the chiral Lebwohl- 
Lasher model, the potential energy is taken as 

= —eP 2 {u> ■ Cl') — K [(u) X (D') ■ r]{uj ■ ui') , (1) 

where ^ 2 ( 2 ;) is the second-order Legendre polynomial. 
The case k = 0 corresponds to the standard, achiral 
Lebwohl-Lasher model. 

These terms are the first scalar and pseudo-scalar 
terms of a general expansion in rotational invariants 
of the interacti on p otential associated with two uniax¬ 
ial particles [1, l36l |. The first, Maier-Saupe-like, term 
favours arrangements where particles are parallel. The 
second term takes account of the chirality of the interac¬ 
tion (in the following the chiral strength will be measured 
by K* = n/e). 

Particles in the Lebwohl-Lasher lattice model can be 
interpreted as single molecules, but Memmer and Flie- 
gans [ 1 ^ have also interpreted each particle on the lat¬ 
tice as a director representing a cluster of neighbouring 
molecules. In our case the different interpretations are 
not relevant considering the qualitative level of our work. 

The setup for the conhned Lebwohl-Lasher model is 
shown in Fig. [1] the lattice is sandwiched between two 
parallel plates normal to the [010] direction. The z axis is 
chosen parallel to the plates, along one of the two equiv¬ 
alent crystalline axes. In the spirit of mean-field theory, 
all particles pertaining to the same plane parallel to the 
surfaces (labelled with the index z = 1, 2, • • • , iV; see Fig. 
[T|) will be taken as equivalent. Particles in planes 1 and N 
interact with their closest plate via the surface potentials 
$i(d>) and $ 2 ( 0 ;). The following cases are considered: 


• Symmetric plates. Both plates have one easy 
axis, contained in the plane of the plates, and the 
two easy axes are identical, n = rh = (0,0,1) 
(see left panel of Fig. [T] for our choice of ref¬ 
erence axes). The corresponding surface poten¬ 
tial energies are monostable, and give contributions 
d>i(cu) = esP 2 {n ■ d>), 1 = 1 and 2, only from the 
particles adjacent to the two plates. 

• Asymmetric plates. The easy axes of the two 
plates form an angle of a = 90°, i.e. n = (0, 0,1) 
and m = (1,0, 0). Contributions from particles ad¬ 
jacent to the plates are $i(tu) = esP 2 {n ■ w) and 
$ 2 ( 0 )) = esP 2 {'fn • tu), respectively. Note that both 
potential energies are monostable, and that their 
strengths are identical; only the easy axis is differ¬ 
ent. 

All surface potentials used are built from second-order 
Legendre polynomials. This is an advantage since then 
all interactions turn out to be quadratic in the unit vec¬ 
tor u), which simplifies the theoretical treatment. We 
take Eg = e (we later argue that this choice leads to strong 
anchoring conditions). In the following we use the Maier- 
Saupe parameter e as an energy scale. We therefore have 
the two parameters kT/e and k/e as energy parameters, 
the length parameter being h/p (both h and p are mea¬ 
sured in lattice unit lengths a, and this is taken as the 
length scale in the model). 

Due to the inhomogeneity of the problem, the orien¬ 
tational order will vary from plane to plane. Since in¬ 
plane fluctuations are neglected, the orientational distri¬ 
bution function will only depend on the plane index z, 
i.e. /i(w), with z = I,-- - , The statistical mechan¬ 
ics of the model is solved using a variational principle 
on the mean-field free-energy functional (for details see 
Marguta et al. El). The functional, has contribu¬ 

tions from entropy, from interactions between particles, 
and from interactions between particles and the plates. 
The energy has contributions from the four particles in 
the same layer, from the two particles on the two neigh¬ 
bouring layers, and also from the external potential. The 
complete free energy functional is: 


Fm] 

eM 


kT J f 

^ N . . 4 

- ^ / dw / dd;7»(‘^)/i(‘^') X! {^2 (w • (1)') + ^ [(‘^ X ‘^0 ■ 7 ] ■ ‘^ 0 } 

i=l j=l 

N-1 „ „ 

^ J du} J dw7*(‘^)/i+i(‘^') {^2 (w • tu') + ^ [(cl) X cl)') • y] (cl) • c2)')| 


+ / dc2)/i(cl))<l)i(c2))-I- / 


( 2 ) 
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FIG. 2. (a) Uniaxial order parameter 77 and biaxial order parameter a (inset) as a function of the reduced temperature fcT/e for 
ft/e = 0.5. The vertical dotted line indicates the isotropic-cholesteric first order phase transition. The curves are qualitatively 
the same for other values of k. (b) Bulk phase diagram in the reduced temperature-K plane. The black-solid line is the 
isotropic-cholesteric binodal. The nematic state is only stable in the limit k ^ 0. The inset is a zoom of the small k region, 
showing the (unstable) isotropic-nematic transition (red line), which does not depend on k. (c) Uniaxial order parameter 77 
and biaxial order parameter a (inset) as a function of k at the isotropic-cholesteric phase transition, (d) Pitch p as a function 
of the chirality parameter k at the isotropic-cholesteric phase transition. The symbols are results obtained via minimization of 
the free energy and the line is the approximation p = 47 r(arctan If )~^. 


where M is the number of particles in each plane, and 
{fj} are the four vectors (±1,0,0) and (0,0, ±1) that 
connect a given particle to its four neighbours in the xz 
plane, and y = (0,1,0) the unit vector along the y axis 
(see Fig. [T]). 

The orientational distribution functions were ex¬ 

panded in six basis functions 'l'a(d)), a = 0, • • • ,5: 

/.(u7) = ^/i“)4/„(±). (3) 

a—0 

The functions 'l'a(w) correspond to real combinations 
of spherical harmonics with angular-momentum indices 
I = 0 and 2 and, therefore, are orthonormal quadratic 
functions in u) (a = 0 corresponds to I = 0 while 
a = 1, • • • , 5 correspond to combinations of the five 


spherical harmonics in the subspace 1 = 2): 

«'o(w) = (4) 

v47r 



4 ' 4 (a>) = Y sin^ 9 cos 2(p, 


4'5 (oj) = Y sin^ 9 sin 2y}. 

The fact that both particle-particle © and particle-plate 
interactions are quadratic in uj limits the (in principle) in¬ 
finite series expansion of fi(pj) to six terms, cf. Eqn. m 
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(the I = 1 subspace is not included due to the head-tail 
symmetry of the particles -no terms linear in u: in the po¬ 
tential). This means that the expansion ([3|) is complete 
and therefore it represents a full, parameterisation-free 
description of orientational order. Since /i(d)) is nor¬ 
malised to unity, = (47r)“^. The five remaining co¬ 
efficients, 

l = (5) 

are a measure of orientational order in each plane i in 
the lab frame. In fact, due to the symmetries of the con¬ 
fined system (the lowest-symmetry phase of which is the 
biaxial phase, with a distribution function fiiCj) which is 
an even function of the angle Lp), the coefficients and 

are zero, and consequently only three order param¬ 
eters, and are relevant. Using Eqn. ([3]) 

the free-energy functional ([2|) can be expressed solely as 
a function of the set of ?>N order parameters for the N 
planes, (in order for this to possible in 

the case of the first, entropic term, the functions 
should be replaced by the expressions following from 
the equilibrium conditions SF/6fi = 0). The necessary 
angular functions were computed using Gauss-Legendre 
quadratures. Minimisation of the free energy was per¬ 
forming using the conjugate-gradient method. The sim¬ 
ple iterative (Picard) method on the equilibrium condi¬ 
tions SF/Sfi = 0 was also used in some cases to check 
the results. 

The order parameters can most easily be interpreted 
in the proper reference frame x'y'z', the one attached to 
the director. Because of symmetry, the director z' lies 
always on the xz plane, at an angle '0i with respect to 
the z axis. The rotation from the lab to the proper frame 
simply involves a rotation about the y axis of angle "0^. 
From the rotational properties of spherical harmonics one 
obtains the relations between the lab, and proper, 

order parameters. In the proper frame only the 

and order parameters are nonzero: 

= P2(cos ^ sin2V’i/f ^ -h ^ sin^ \ 

^(4) ^ ^ 1 ) -h P2(cosV'i)/f’(6) 

is related with the standard uniaxial order parame¬ 
ter, rji = which measures the uniaxial order 

of particles about the local director. The other param¬ 
eter is related with the usual biaxial order parameter, 
ai = Y^7r/15, which accounts for possible departures 

from uniaxial symmetry of particles about the local di¬ 
rector. The third and last relevant parameter is the tilt 
angle tiji (which also depends on the plane index, since in 
the cholesteric phase the director rotates from plane to 
plane). 



FIG. 3. Uniaxial (a) and biaxial order parameter (b) profiles 
of an isotropic (open squares) and a cholesteric (filled circles) 
phase coexisting in a symmetric slit pore with N = 25. The 
tilt angle of the cholesteric phase (not shown) rotates by ap¬ 
proximately 27r across the cell. 

III. RESULTS 
A. Bulk system 

We first describe the bulk phase behaviour of the 
model. In the achiral system (k = 0) there is an isotropic- 
nematic phase transition (s^, l3^ . For any nonvanish¬ 
ing chirality {k ^ 0) the nematic phase is replaced by 
a cholesteric phase, and the system exhibits a first or¬ 
der isotropic-cholesteric (Iso-Cho) phase transition by 
decreasing the temperature. The uniaxial and biaxial 
order parameters as a function of the temperature for 
At/e = 0.5 are shown in Fig. Ha). The cholesteric pitch 
(not shown) is inde pen dent of the temperature, as can 
be demonstrated [^, [ij . The behaviour of the order pa¬ 
rameters is qualitatively the same for any value of k. 

In Fig. Ht>) we represent the bulk phase diagram in 
the plane of chirality and temperature. The stronger the 
chirality is, the higher the phase transition temperature 
becomes. In the limit of vanishing chirality (k —a 0), 
the transition temperature tends asymptotically to the 
isotropic-nematic transition temperature of the achiral 
system. 
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The uniaxial and biaxial order parameters at the phase 
transition are depicted in Fig. He). As expected the 
biaxiality is very small. Finally, in Fig. we show 
the cholesteric pitch at the Iso-Cho phase transition as 
a function of k. The cholesteric pitch increases rapidly 
for decreasing chirality, diverging (p —^ oo) in the limit 
K —?> 0. This limit corresponds to a system of achiral par¬ 
ticles in which a nematic or, equivalently, a cholesteric 
with infinite pitch, is stable. The data closely follow 
the curve p = 47r(arctanan approximation that 
results by neglecting the entropy and minimising the in¬ 
ternal energy with respect to the pitch. 


B. Confined system 

Next we confine the system in a symmetric pore made 
of two identical parallel plates. As mentioned before, we 
set the surface strength as Cg = e. This choice corre¬ 
sponds to moderate or strong anchoring conditions on 
the two plates since, for relatively thin pores, the parti¬ 
cle orientation next to the plates never exceeds 2° from 
the easy axis for stable configurations. Although we have 
studied different values of k, here we only show results 
for k/c = 1, as no qualitative differences were found for 
different values of the chirality strength. For this choice 
of K the cholesteric pitch is p ~ 21 (in units of the lattice 
parameter). 

1. Symmetric cell 

In the symmetric cell the anchoring directions are 
the same, i.e. n = m (both vectors lying on the 
plane of the plates). In pores narrower than approxi¬ 
mately half a cholesteric period {N < p/2), the capil¬ 
lary isotropic-cholesteric transition disappears, and the 
cholesteric grows from the isotropic state in a continu¬ 
ous fashion as temperature is decreased. The suppres¬ 
sion of the capillary transition for narrow pores is typical 
of confined fluids. In a narrow pore surface interactions 
strongly affect the particle orientations in all layers, not 
just in those adjacent to the plates. Bulk orientational 
properties cannot compete with the strong anchoring im¬ 
posed by the plates, and the absence of the bulk-surface 
competition results in a perfectly continuous evolution of 
orientational order parameters, as dictated by surface po¬ 
tentials and temperature. Despite the fact that there is 
no Iso-Cho phase transition, the system director, whose 
(symmetric) profile is close to linear, does show a total 
twist mainly dictated by pore width and nearly indepen¬ 
dent of temperature. The total twist decreases as the 
pore becomes narrower. 

In wider pores {N > p/2) we found a first-order Iso- 
Cho phase transition. Order parameter profiles of the 
confined isotropic and cholesteric phases coexisting in a 
pore with N = 25 are plotted in Fig. [31 In the isotropic 
state the uniaxial order parameter, panel (a), vanishes in 


N/p 



N/p 


FIG. 4. Phase diagram of a confined cholesteric liquid crystal 
in (a) symmetric, and (b) asymmetric pores in the scaled- 
temperature T/T* versus scaled pore-width A/p plane. T* is 
the temperature for the bulk isotropic-cholesteric phase tran¬ 
sition. Filled circles indicate the isotropic-cholesteric bin- 
odal. Dashed lines represent the binodals along which two 
cholesteric phases with different number of helix indices co¬ 
exist. The dotted red line indicates the temperature of the 
bulk transition. In each case the geometry of the easy axes is 
depicted in the upper-right corner of the corresponding panel. 

the middle of the pore, and it is relatively high close to 
the plates. In the cholesteric phase, the uniaxial order 
parameter is always higher than zero, and particles are 
also more ordered close to the plates. The biaxial order 
parameter, panel (b), is very small in the whole cavity. 
As can be seen, it is always negative, meaning that parti¬ 
cles fluctuate about the director with a larger amplitude 
in the direction of the helix (plate normal) than in the 
direction of the plates. 

The temperature at which the Iso-Cho transition oc¬ 
curs varies with pore size A. This is reflected in the phase 
diagram, depicted in Fig. jH)^) ^ ^^e pore size (scaled 
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FIG. 5. Total twist tp (in degrees) of three cholesteric phases 
in a symmetric pore as a function of pore width scaled with 
cholesteric pitch, N/p. From bottom to top: Cstt and 
phases. In each case results are given for three different tem¬ 
peratures: continuous lines T/T* = 0.800 (where T* is the 
bulk transition temperature), dashed lines TjT* = 0.880, and 
dotted lines T/T* — 0.987. Vertical dotted lines correspond 
to the C^,TT — CeTT phase transition, whereas dots indicate the 
corresponding coexisting phases. 


with the cholesteric pitch) versus temperature (scaled 
with the bulk Iso-Cho transition temperature) plane. We 
first focus on the capillary Iso-Cho transition line. The 
transition is shifted to higher temperatures with respect 
to the bulk; therefore, conhnement promotes cholesteric 
order. The capillary line oscillates and tends asymptoti¬ 
cally to the bulk transition temperature. The oscillatory 
behaviour is a direct consequence of a commensuration 
effect between the pore size and the cholesteric period. 
Local maxima of the capillary line occur when the elastic 
stress of the confined cholesteric is minimal, which hap¬ 
pens if the pore size and the cholesteric period commen¬ 
surate, i.e., N/{p/2) ~ k, with k a positive integer (note 
that, due to the head-tail particle symmetry, the period 
of the cholesteric phase is in fact p/2, corresponding to 
half a turn of the director about the helix axis). Any¬ 
way the scaled-temperature range where this oscillatory 
structures takes place is rather small. 

Besides the capillary Iso-Cho line, the system exhibits 
winding transitions in which two cholesterics with differ¬ 
ent periods coexist. The winding transitions occur always 
between cholesterics that differ in half a cholesteric pe¬ 
riod: CiTT — where the subscript in, i = 1,2,3..., 

denotes the total twist of the cholesteric helix from one 
plate to the other (save the small deviations of the di¬ 
rector at the two plates arising from a large, but fi¬ 


nite, anchoring strength). At a winding transition, the 
cholesteric period of the coexisting phase is expanded 
with respect to its equilibrium value, while the period 
of the coexisting phase is compressed. Hence, 

the size of the pore at which two cholesterics coexist is 
such that the elastic stress is maximum. Winding tran¬ 
sitions merge with the Iso-Cho capillary line, giving rise 
to Iso—Ci^ — triple points. The triple points are 

located on the local minima (low temperature) of the 
capillary Iso-Cho line. 

Note that winding transitions are almost vertical lines 
in the pore size-temperature phase diagram, implying 
that they are insensitive to temperature. This result can 
be traced back to the insensitivity of the cholesteric pitch 
with respect to temperature. To see this more clearly. 
Fig. Oshows the total twist across the cell (from one plate 
to the other) of the director as a function of pore width 
N/p for two phases, and Cqt^, in the region of the 
Cstt ~ Cqtt phase transition. For each phase, three curves 
corresponding to three temperatures, T/T* = 0.800, 
0.880 and 0.987 are shown. The three curves (and, in 
fact all curves corresponding to different temperatures) 
intersect at a point with a total twist equal to the num¬ 
ber of half-periods, Stt = 900° for the phase and 
6n = 1080° for the Cqtt phase, and at a pore width 
slightly higher than the optimal value (which is given by 
N = 5{p/2) = 2.5p and N = 6{p/2) = 3p, respectively); 
this is due to the effect of the plates, which disappears as 
the pore width becomes larger. The — Cstt transition 
is obtained when the free energies of the two structures 
become equal. For all three temperatures, the transition 
is located at almost exactly the same pore width (denoted 
by the dotted vertical line in the figure). In effect, due 
to the temperature invariance of the pitch, the tempera¬ 
ture plays a negligible role in driving winding transitions. 
By contrast, a model with temperature-dependent pitch 
would give rise to boundary lines presenting a slope in 
the phase diagram. Finally, as mentioned before, we note 
that, at coexistence, the phase with the smaller number 
of turns, Cstt) is expanded, whereas the one with the 
larger number of turns, Cqtt, is compressed. This is be¬ 
cause the transition point is located at wider (narrower) 
pores than optimal for the Cs^r (Ceir) phase. 

2. Asymmetric cell 

Next we confine the system in an asymmetric pore, 
in which plates favour perpendicular directions (easy 
axes satisfying n T m). All remaining parameters 
are set to the same values as in the previous symmet¬ 
ric case. Not surprisingly the phase diagram, shown 
in Fig. SKb), resembles the one of a symmetric pore. 
The only significant difference concerns the total twist 
of the confined cholesterics, which is always a multiple of 
7r/2; therefore the cholesteric phases can be labelled as 
C( 2 j-i-i)ii -/2 = C'(i+i) 7 r) with i = 0,1,2.... The subscript 
{i + ^)n indicates the total twist inside of the pore. Again 




cholesteric phases differing in half a cholesteric period, 
~ coexist at the winding transitions. 

The phase diagram is essentially shifted by half a period 
with respect to the symmetric case. However, the pore 
width below which the Iso-Cho capillary line disappears 
is very similar to that in the symmetric cell. 

IV. CONCLUSIONS 

In this work we have calculated the phase diagram of a 
cholesteric liquid crystal confined in a slit pore between 
two parallel plates. To our knowledge this is the first time 
the structure of winding transitions is explicitly obtained 
as a function of pore size and connected to the capillary 
transition. We have investigated symmetric (identical 
plates) and asymmetric pores (easy axes of plates at an 
angle) pores. As a result of confinement, the cholesteric 
phase is promoted with respect to the bulk, the binodal 
line exhibiting a oscillatory structure, and strong coin- 
mensuration effects, giving rise to winding transitions, 
occur. The two types (capillary and winding) of transi¬ 
tions are connected: the capillary line preempts the wind¬ 
ing transition structure, and the oscillatory nature of the 
former is due to the commensuration effects at work in 
the winding transitions. Our results pertain to the case 
where the pore size of the cell is changed in the absence 
of an external field, a situation not common in practical 
applications; however, these effects have been observed in 
Atomic-Force Microscopy experiments, where both pore 
size and temperature can be controlled. 


Finally, it is interesting to note that phase diagrams 
similar to the ones obtained here have been previously 
predicted in smectic (iol . 1^ or columnar [42j liquid crys¬ 
tals confined in slit pores using density functional the¬ 
ory. In these system layering transitions between states 
that differ in the number of smectic or columnar layers 
occur under confinement due to commensuration effects 
between the pore width and the layer spacing. Layer¬ 
ing transitions play a role similar to the winding transi¬ 
tions analysed here; they are connected to the capillary 
isotropic-smectic or nematic-smectic lines, a phenomenon 
which is the analogue of the connection between the 
winding transitions and the capillary isotropic-cholesteric 
line found for the cholesteric model. This analogy results 
in a universal behaviour, stemming from the competi¬ 
tion between an internal length and an imposed, exter¬ 
nal length. However, some differences exist. For exam¬ 
ple, the capillary transition lines at small pore widths is 
in some cases broken in the smectic case. Possible ori¬ 
gins of these differences (resulting from the theoretical 
approximations or from the interaction models) is cur¬ 
rently under investigation. 
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